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We have investigated the relation between the two-point spatial correlation 

function and the density profile of the dark matter halo in the strongly non-linear 
regime. It is well known that when the density fluctuation grows into dark 
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matter halo whose density profile is p oc r e (| < e < 3) on almost all mass 



scales, the two-point spatial correlation function obeys a power law with the 
power index 7 = 2e — 3 in the strongly non-linear regime. We find the form of 
the two-point spatial correlation function, which does not obey the power law 
when the power index e is smaller than |, such as the density profile p oc r _1 
Q\ . around the center of the halo which is proposed by Navarro, Frenk & White 

(1996,1997). By using the BBGKY equation in the strongly non-linear regime, 
it is also found that velocity parameter h = — {v)/ax is not a constant even in 
the strongly non-linear regime (x = x/x n i — > 0) although it is a constant when 
e > 3/2 and then the two-point spatial correlation function can be regarded as 
the power law. The velocity parameter h becomes at the non-linear limit of 



Or 



x — > 0, that is, the stable clustering hypothesis cannot be satisfied when e < 3/2. 
Subject headings: cosmology:theory-large scale structures-dark matter halo 

1. INTRODUCTION 

Formations of the large scale structures in the expanding universe is one of the 
most important and interesting problems of cosmology. It is generally believed that 
these structures have been formed owing to the gravitational instability. Hence it is very 
important to clarify evolutions of density fluctuations by the gravitational instability. Here 
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we consider density fluctuations of collisionless cold particles such as cold dark matters, our 
interest being mainly concentrated on effects of the self-gravity 

In the hierarchical clustering picture, small fluctuations grow and form small clusters 
of dark matter (dark halo) at first. Then these small dark halos merge into a larger halo 
or accrete surrounding matters and produce the larger dark matter halos. Evolutions of 
density fluctuations in the non-linear regime are especially important to understand the 
formation of the dark matter halos. 

A central region of the dark halo profile was investigated by using N-body simulations 
in 1990's (Frenk et al. 1988, Navarro, Frenk & White 1996, 1997 : hereafter NFW ). 
Hernquist (1990) investigated the elliptical galaxy analytically. NFW and Hernquist 
claimed that dark halos are not well approximated by isothermal spheres. Instead, the 
density profile is well approximated by the following form: 



p ° Vhi + ^ ' (1) 

v r s ' v r s ' 

where, r s is a characteristic scale of the dark halo. The parameter, e, and /i are equal to 1 
and 4, respectively, in the Hernquist's case, and e = 1, /i = 3 in the NFWs case. In any 
case, the density profile of the central region is shallower than the isothermal sphere. 

On the other hand, the two-point spatial correlation function is well used in the analysis 
of the density fluctuations of the dark matters. Of course, the two-point spatial correlation 
function in the non-linear regime is related to the dark halo profile if the mass weight of the 
dark halos dominates in the universe, that is, almost all regions of high density correspond 
to dark halos. Davis & Peebles (1977) assumed the self-similar evolution and the stability 
condition, that is, the mean relative physical velocity (r) is equal to in the non-linear 
regime. Then, they showed the index of the two-point spatial correlation function in the 
non-linear regime by using the index of the initial power spectrum, n: 

c -7 3(3 + n) 

£ocr 7 , 7 = — . (2) 

5 + n 

When we do not assume the stability condition, the index of the two-point spatial 
correlation function 7 is expressed by using the relative velocity parameter h = —{v)/ax, 

3M3 + n) ,o\ 

7 2 + h(3 + nY [) 

where n > —3, and h is a value of order 1 (0 < h < 1) (Yano & Gouda (1997) ; hereafter 
YG) in the non-linear regime. Here (v) is the mean relative peculiar velocity ( see eg. (|T0D ) 
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and a is the scale factor. Therefore 7 have the value around the following: 

< 7 < 3. (4) 

We consider the case that the halo dominates in the regions of high density and 
then the density profile of the halo mainly contributes to the two-point spatial correlation 
function in the non- linear regime. When the density profile of the dark halo is given by 
p oc r~ e , then using the following equation, 

£(r)~(p( S )p(|r + l?|)>/<p> 2 , (5) 

the two-point spatial correlation function in the strongly non-linear regime (r — ► 0) is 
determined by the density profile of the dark halo and written in the following way 
(McClelland & Silk 1977; Sheth & Jain 1997, Padmanabhan & Engineer 1998), 

£ocr" 7 ', y = 2e-3. (6) 

The two-point spatial correlation function diverges at the center of the dark halo when 
e > 3 and it does so on r = r/r n i — > 00 when e < 3/2, where r n \ is the non-linear scale 
of the density fluctuation. Therefore 3/2 < e < 3 must be satisfied. When the relation 
3/2 < e < 3 is satisfied, 7' satisfies < 7' < 3. This relation corresponds to the relation (||) 
derived from eq. ([3]). In other words, when the two-point spatial correlation function obeys 
the power law with < 7 < 3, the index e satisfies 3/2 < e < 3. 

When the density profile is p oc r~ e with e < 3/2 around the center of the dark halo, 
the index of the density profile have at least following two values in order to converge the 
integral. 

p oc (L)-M : <r) = J e <T < '. : e < I) (7) 

r/ { p (r>r s : // > §) 

Here r s is a characteristics scale. In this case, what is the form of the two-point 
spatial correlation function ? The purpose of this paper is to obtain the two-point spatial 
correlation function in the strongly non-linear regime when the dark halos have the above 
density profile (0) and the self-similar evolution of the dark halos is satisfied. We also 
investigate what is the form of the velocity parameter h, that is, the mean relative peculiar 
velocity by using the BBGKY equation. 

In §2, we will briefly review the power law solution of the two-point spatial correlation 
function described in YG. We show in §3 the solution of the two-point spatial correlation 
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function and the mean relative peculiar velocity when there exists the dark halo whose 
density profile has the form of eq. (|7|) . We will devote to the conclusions and discussions 
in §4. In the Appendix, we will consider the relation between the power index of the 
density profile of the dark halo and the initial power index n if the merging or accreting far 
equilibrium dark halos could be regarded as a spherical halos. 



2. POWER LAW SOLUTIONS 

We will briefly review the power-law solutions of the two-point spatial correlation 
function £, and the mean peculiar velocity (v) described in YG. We use the second BBGKY 
zeroth moment equation. 

where, a is the scale factor, x is the comoving coordinate. The mean relative peculiar 
velocity vector (v) is defined as follows; 

, , = (S(V2 ~ v 1 )f{l)f{2)dh 1 dS 2 ) 

{ } ~ (J f{l)f{2)dH x d?v 2 ) • l ] 

Here (...) shows the ensemble mean for any pair with the fixed distance x = \x 2 — X\\ 
and f(i) = f(xi,Vi) (i = 1,2) is the phase space density at the point of the phase space 
(xi,Vi). We notice that (/(l)) = 6(1) is the one-body probability distribution function, 
(/(l)/(2)) = 6(1)6(2) + c(l, 2) is the two-body probability distribution function, and c(l, 2) 
is the irreducible two-point correlation function. Here we assume that the background 
universe is homogeneous and isotropic. Then, (v) oc x 2 — x± = x and then the scalar 
velocity (v) = (v)— is the only important quantity which determines £ through eq.(||). The 
mean relative peculiar velocity (v) can be rewritten by 

Therefore we can define the mean relative peculiar velocity (v) as follows: 

(f_f ■(v 2 -v 1 )f(l)f(2)d 3 v 1 d^ 2 ) 
W Uf(l)f(2)dz Vl d3v 2 ) ■ l } 

We assume that the two-point spatial correlation function £ is given by the following 
power-law form in the strongly non-linear regime: 



^ = ^x~\ (11) 
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Then we obtain from the dimensional analysis of the equation (||) 

(v) = -hax, /3 = (3-7)/i, (12) 

where h is a constant. We call this parameter h, the velocity parameter. This parameter is 
a value of order 1 (0 < h < 1)(YG) when we consider the collisionless cold dark matters in 
the strongly non-linear regime. We notice that the stability condition proposed by Davis & 
Peebles (1977) corresponds to h — 1 (stable clustering). In this case, the collapsed object 
cannot be broken and clustered together to form the larger cluster. Therefore, the relative 
physical velocity of two particles (f) = (v) + ax is equal to ((v ) = —ax) in the non-linear 
regime. We can consider the other extreme case about the clustering picture. In this case, 
the smaller objects have clustered and merged together and the completely virialized object 
is newly formed. In this case, the separation of the particle is expanding with the Hubble 
velocity and then h = 0(comoving clustering). 

If the self-similarity solutions exist, the power index of the two-point spatial correlation 
function in the non-linear regime can be represented by the following form (Padmanabhan 
1996; YG): 

3h(3 + to) 

2 + h(3 + to) 

where to is the index of the initial power spectrum (n > —3). Therefore 7 have the value 
between and 3. Here we consider the case that the halo mass dominates and the index 
of the density profile p oc r~ € is the same for all mass scales. The relation 3 — e > and 
3 — 2e < must be satisfied so that the integral in eq.(]5|) should converge. Therefore, e have 
the following values: 

\ < e < 3- (14) 

In this case, the relation between the two-point spatial correlation function and the density 
profile is 



£(r) oc / d 3 sp(s)p(\l^ + if |) 

= c 2 (— )" 2e+3 . (15) 

r n i 

[see McClelland & Silk 1977; Sheth & Jain 1997, Padmanabhan & Engineer 1998] where, 

00 1 1 1 

c 2 = V 2( )(2 - e)(l - e) • • • (2 - e - 2m)- -. (16) 

^ V 3-e + 2m 3-2e-2m A A J K ; (2m+l)! y J 



-6- 



The relation between the both power indexes of the halo profile e and the two-point 
spatial correlation function 7 is given by 



7 = 2e-3. (17) 

We can see from equations (p~^|) and (|T7j ) that the index of the two-point spatial correlation 
function 7 have the values between and 3. This result is consistent with one derived from 
eq.©. 



On the other hand, when we assume the density profile p oc r~ e with e < 3/2 around 
the center of the dark halo, what is the form of the two-point spatial correlation function? 
We consider the relation between the two-point spatial correlation function and this density 
profile in the next section. 



3. UNIVERSAL DENSITY PROFILE 

We consider the form of the two-point spatial correlation function when the density 
fluctuations grow into the dark halos which have the shallower density profile p oc r~ e with 
e < 3/2 around the center. In this case, the density profile cannot have the single power 
law with e < 3/2 in order to converge the integral about the two-point spatial correlation 
function in the r = rjr n \ — > 00 limit. 

Navarro, Frenk & White (1996,1997) and Hernquist (1990) claimed that the halo 
density profiles have the following form: 



P0C (^(1 + ^ ' (18) 

v r s ' v r s > 

where p > 3/2 must be satisfied because the integral of the density should converge (eq. 
(|l8l)). The characteristic scale r s can be rewritten by Xr n i, where A is a constant. For easy 
treatment, we assume that the density profile is written in the following way: 

K(r) . .^ _ j e (r < Xr nl : e < 3/2) 



'^-W"" : «M= „ r>*: ;„>3/2 ■ < i9 > 



Here we consider the case of A < 1 and normalize the density such as p = 1 at r = r n i 
although the same results are obtained in general for A > 1. The two-point spatial 
correlation function is calculated by the following form: 
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£(r) oc / d 3 sp(s)p(|T> + ^|). (20) 



Here we have assumed that the halos are formed in the non-linear regime and their 
density profiles have the same form independently of mass scales as we show in eq. (|H?D , 
that is, we have assumed that their halos evolve self-similarly because many authors such 
as Jain & Bertschinger (1996), Colombi, Bouchet, & Hernquist (1996), Yano & Gouda 
(1998) have investigated the self-similarity about the two-point spatial correlation function 
and then the self-similarity is widely believed. In this case, each halo should have the 
self-similarity in the same way as the two-point spatial correlation function in the non-linear 
regime. Furthermore, the correlations with the other halos do not affect the value of the 
two-point spatial correlation function in the non-linear regime because we consider the 
non- linear regime (x = x/x n i — > 0, where x n i is the non-linear scale written in the comoving 
coordinate), and then almost all pairs of the two points with the distance x <C 1 are 
included in one halo. Therefore, it is enough to consider the two-point spatial correlation 
function within only one halo in order to obtain the correlation function at x <C 1. Then, 
the two-point spatial correlation function in the non-linear regime can be derived as follows: 



f(r) oc / d 3 sp(s)p(\ ~r> + -?\) 



OC ^1 [ dsdcos s 2(_!_)-«(r)(JL)-*(v) 

Ki Jhaio \r nl Xr n i 

- c 2 (-)- 2e+3 + tt^A 3 "^ + ^A 3 - 2 - + 0((^) 2 ) 
r ni 3 - 2e 2/i - 3 r nl 

= c 1 + c 2 (r) 3 - 2e + 0(r 2 ) (21) 



where y = \/r 2 + s 2 + 2rs cos 6, 9 is the angle between the direction of "r* and if, c\ and 
C2 are constants, and c\ = A 3 ~ 2e /(3 — 2e) + A 3_2M /(2yU — 3) is a positive value because we 
consider the case that e < 3/2, and /x > 3/2. As a result, we can see that in the case 
of e < 3/2, the two-point spatial correlation function cannot obey the power-law in the 



non-linear regime. When the density is continuous at r = r s as shown in eq.(|l^), the term 
of 0(r) is eliminated. On the other hand, the term of 0(r 2 ) remain. However, there is 
the possibility that the term of 0(f 2 ) or the higher terms are eliminated when the density 
profiles are smoothly continuous at r = r s . In the following, we investigate the parameter h 
when the two-point spatial correlation function can be expressed by eq.(pTl). 



Here we put the two-point spatial correlation function £ and the velocity parameter h 
in the following way: 



£ = C\ + c 2 a 7a x 7 

/i = t/i + ^a-^F, (22) 

where x = x/x n i, and a is a constant value for the scaling. Here ci, C2, di, and <i 2 are 
constants. When the self-similarity of the two-point spatial correlation function is satisfied, 
a is equal to 2/(n + 3) (DP, YG). The index 7 satisfies —7 = 3 — 2e, (e < 3/2) (see eq.(PT])). 
We define the velocity parameter in the same way as the power-law case ((v) = —hax). 
Then the BBGKY equation (§) reduces to 

%ac 2 a ai x^ = - — — [ x 3 ( Cl + c 2 a fa x~ f ){d 1 + d 2 a~' MX x 1 ")], (23) 



a 



ax 2 dx 



and then, this is rewritten as follows: 



7 ac 2 a Q7 5- 7 = 3adi + {v + S)c x d 2 a- V<x x v + (-7 + 3)dic 2 a 7Q a;- 7 

+ (-7 + z/ + 3)c 2 d 2 a( 7 ~^s- 7+ ^. (24) 

Since we consider e < |, the relation —7 = 3 — 2e > is satisfied. Therefore, the left 
side of equation (^|) is in the non-linear limit [x = x/x n i — > 0). As a result, d\ must 
be because c\ is a positive value. Then the first and third terms of the right side of the 
equation (|4]) vanish. Therefore in the limit of x — > 0, the second term of the right side is 
the same order of the left side. So, v is equal to —7. Then, the last term of the right side is 
higher order. Thus we obtain the following relation 



7«c 2 = [y + 3)ci<i 2 , (25) 

or equivalently, 

-3cid 2 , „x 

c\d 2 + ac 2 

Here a and c 2 are order of 1. 7 and v + 3 are the same order. c\ is the value of the 
two-point spatial correlation function in the non-linear limit (£(0)). Then c\ is greater than 
1. If we assume C\ ~ 10 2 , then d 2 ~ 10 -2 . 

Finally, we find that the velocity parameter is written in the following form: 

h = d 2 ar va x v : z/ = 3-2e (e < -), d 2 < 1. (27) 
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The velocity parameter h is zero in the non-linear limit. This means that the stability 
condition (h = 1) cannot be satisfied at least in the case of e < 3/2. 



We obtain the form of the velocity parameter as shown in eq. (j27|) , and this velocity 
parameter is related with the two-point spatial correlation function through v. However, we 
can not obtain the values of d 2 and v from the second BBGKY zeroth moment equation. If 
we would like to do so, we must solve higher moment equations. 

Here we consider the meaning of the velocity parameter h again. The parameter h is 
defined as follows: 

W = -*" = (//(1)/(2)A 1 A,> • (28) 

This velocity parameter is derived from an ensemble mean of the relative peculiar velocity, 
and can be regarded as a spatial mean of all regions of the universe when the distance of 
the two-points is much smaller than the horizon scale and we can get many samples of the 
pairs. Then if we average the relative peculiar velocity within a part of the region of the 
universe, this value is generally different from the mean value averaged within all regions of 
the universe. Here we define the following parameter h\ oca i\ 

{v) local - h local ax- (Jf(l)f(2)d^v^v 2 ) locai ' [ZJ) 

where the subscript "local" means that the average is taken only in a dark halo for any 
pair with the fixed distance x = \x 2 — x±\. When a halo becomes isolated and virialized, 
this halo might have the value of the parameter hi oca i = 1 even when the stability condition 
cannot be satisfied (h ^ 1) (see Appendix). It should be noted that the parameter hi oca i = 1 
does not necessarily mean h = 1. 

We obtain the result of h = in the non-linear limit as we can see from eq. (ffi\)- Then 
if hi oca i > in a dark halo, there must exist another halo in which hi oca i < 0. Because if all 
halos have hi oca i > 0, h should be larger than and it is inconsistent with the condition of 
h = 0. 

Here we assume that the density profile is determined independently of the initial 
conditions. In this case, a kind of "relaxation" process must work for the density profile. 
Then, the velocity field of each halo also might relax to a common state independently of 
the initial condition. Therefore, all halos might have the same hi oca i. Since we obtain h = 
in the nonlinear limit, hi oca i for all halos must be 0. 

Only if the power index e of the density profile is equal to 0, the both conditions of 



-10- 



h = and hi oca i = for all halos can be satisfied (see eq.(|39|) in Appendix). If e ^ 0, as 
appeared in the NFW density profile, then hi oca i ^ 0. So it is impossible that all halos with 
e^O have the same value of hi oca i ^ independently of initial conditions. 



4. CONCLUSIONS AND DISCUSSION 

We have investigated the relation between the two-point spatial correlation function 
and the density profile of the dark matter halos in the strongly non-linear regime. It is 
found that when the density fluctuations evolve to the dark matter halos whose density 
profile is p oc r _e (e < 3/2) around the center of the halo, the two-point spatial correlation 
function cannot obey the power law. Furthermore we find by using the BBGKY equation 
that the velocity parameter h = —(v)/ax is not a constant in the case of the shallower 
density profile (e < 3/2) although h is a constant in the case of e > 3/2. The velocity 
parameter h becomes at the non-linear limit, that is, the stable clustering hypothesis 
cannot be satisfied at least in the case of e < 3/2. 

We consider the power spectrum when the index of the density profile is shallower than 

3/2: 



P(k) = U(r)e ikx d 3 r 

ft i r x,,! sin £;r 



r dr 



3 3 + z/ v ; 

The index v = 3 — 2e of the power spectrum is 

dlnP_ vKk~» 

d\nk~ 6 l + Kk-» [> 6h [6l) 

where K = 3c2r^/(3 + u)ci(< 0). In the case of the power index ^^ < —3, the power 
spectrum evolves and the power index changes to the value —3 which is decided by the 
catastrophe theory (Gouda 1998). However, as we see from eg. (pl|) , the index of the power 
law is larger than —3 when the power spectrum is fitted by the single power law in large k 
although the power spectrum does not obey the power law in real. 

In the paper of the NFW, it is proposed that the density profile is independent of the 
initial condition and the index of the density profile e is equal to 1. Here we assume that 
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the NFW density profile is realized in almost all halos on all mass scales of the halo. In this 
case, a kind of "relaxation" process must work for the density profile. Then the velocity 
field also might be relaxed in all halos. Therefore all halos might have the same hi oca \. Since 
we obtain h = in the non-linear limit, the parameter h\ oca i must be if hi oca i is common 
for all halos. As we can see from the argument in Appendix ( eq.(|3~9"l)), hi oca i becomes 
only when the index of the density profile e = 0. Primack & Bullock (1998) proposed by 
simulations that the density profile in the center of the dark halos have the index with 
e ~ 0.2. This value is different from the value proposed by NFW. Furthermore this value 
e p» 0.2 is nearly equal to 0. Therefore there exists the possibility that all halos have the 
same hi oca i and its value is 0. If the index of the density profile of the dark halos have the 
value of e = 1 that is proposed by NFW, hi oca i must not be in the central region of the 
dark halo. In order to become h = 0, there might exist another halo in which hi oca i < 0, 
that is, each halo might have a different velocity field. Here it is hard to believe that each 
halo has independently the same index with e = 1 although each halo has the different 
velocity field. Then, it is natural that the index e also depends on the initial condition or 
e = for all halos on all mass scales. 

On the other hand, Huss, Jain, & Steinmetz (1998) have claimed that the density 
profiles of the isolated object have a nearly universal characteristic shape independently 
of the initial condition or the formation history. This means that the halo which is saved 
from the merging or the accretion for enough time to become equilibrium might have the 
index e = 1. If the dark halos with the index e = 1 whose density profile is determined 
independently of the initial condition dominate in the universe, each halo might have the 
same hi oca i and hi oca i is larger than (see eq.([39|) in Appendix), and then h cannot be equal 
to 0. Therefore, such halos can not dominate because h must be in the non- linear limit. 
In order to become h = 0, we can consider that such virialized and equilibrium halos do not 
dominate and far equilibrium and distorted halos dominate in the universe. And then, the 
far equilibrium halos, that is, the merging or accreting distorted halos mainly determine the 
self-similarity of the two-point spatial correlation function. 

We consider the case that almost all halos are merging or accreting and are not 
equilibrium. Such halos can not be spherical in a real situation. However, we assume the 
distorted dark halos could be regarded as spherical halos. Then we can regard p as a function 
of r only. In this case we investigate the density profile of the dark halos, the two-point 
spatial correlation function, and the velocity parameter as a function of the index of the 
initial power spectrum as we show in the Appendix. We assume that the self-similarity of 
the dark halos satisfies and that hi oca i is equal to 1. In this assumption, the index of the 
density profile becomes e = 3(n + 3)/(n + 5) from eq.pO]) in the Appendix. In the case of 
n > — 1, that is, the index of the density profile is 3/2 < e < 3 because e = 3(n + 3)/(n + 5). 
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Furthermore the velocity parameter is h = (n + l)/2(n + 3) and the index of the two-point 
spatial correlation function is 7 = 3h(n + 3)/[2 + h(n + 3)] = 3(n + l)/(n + 5). In the case 
of n < — 1, the index of the density profile is e < 3/2. In this case, the power index of the 
two-point spatial correlation function in the non-linear limit is (£ = C\ + C2(r/r s ) 3-2e ), 
and the velocity parameter in the non-linear limit is (h — d,2X i ~ 2e ). 



5. APPENDIX 

The merging or accreting dark halos are distorted and far equilibrium. In this 
Appendix, we show the relation between the power index of the density profile of the halo 
and the initial power index n if such distorted dark halos could be regarded as a spherical 
halos and we can obtain the density profile which depends only on the distance r from the 
center of the halo. When we assume that the self-similar evolution and the hierarchical 
clustering are satisfied, the density of the dark halo can be written as follows: 

p(r) = M-r, (32) 

r s 

where 5m oc a~ 3 is the normalization factor and r s oc a 2 ^ n+3 ^a = a^™ +5 - ) ^™ +3 - ) is the 

characteristic size of the dark halo (Peebles (1980), Syer & White (1997)). The continuity 

equation of the matter in the halos is given by 

^ + V • (pV) = 0, (33) 

where V is the physical streaming velocity. 

}Vf(2)dh 2 

where V is the physical velocity and /(2) is the phase space density. Here we express the 
streaming velocity by using a newly defined parameter h stream 



V r = V ■ r = v(x) + ax = (1 - h stream )ax, (35) 

where x is the separation from the center of the dark halo and v(x) is the peculiar streaming 
velocity. Hence the parameter h str eam is expressed by the following form, 

., _ Of ■V2f(2)d 3 V 2 ) x 

an stream x ^f(2)d?v 2 ) x ' l ' 
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where (...) x means the average in a dark halo for the points which have the fixed distance 
x from the center of the dark halo. If the radial streaming velocity around the center 
of the dark halo is vanished due to the virialization, the following relation is satisfied; 
(V r = V ■ r = 0). Then in this case h stream is equal to 1 as we see from eq. 



The peculiar streaming velocity v corresponds to the mean peculiar velocity of the two 
points, one of which is fixed to the center of the dark halo. Therefore v is different from 
(v) local, which is defined in eq. (p9|) , in general. However, there is the possibility that v is 
equal to (v)i oca i, that is, h stream is equal to hi oca i- For example, we consider the case of 
Vi oc Xi where Xi is the radial vector of the i-th point from the center of the halo. This 
relation might be satisfied when the halo is isolated and virialized. In this case we can show 

hstream = hocal aS follows: 



ahi oca ix 



(I f • (V2 ~ Vl)f(l)f(2)^Vl<Pv 2 )local 

(J f(l)f(2)d 3 Vl d s v 2 ) local 

(- • (v 2 - Vl)) local 

(- ■ (v a -vi))x 

X 

(- ■ v 2 ) x - {- ■ v 1 ) x 

(- ■ v 2 ) x 
(Jf-v 2 f(2)d 3 v 2 ) x 



(ff(2)d*v 2 ) a 



-ah 



stream 



X 



(37) 



In the above equations we use the following arguments; Since Vi oc x^ the relation 
v 2 — Vi oc x 2 — x\ = x is also satisfied for any two-points in a halo. Therefore, if we 
consider the average for the fixed x±, then the mean value of the relative peculiar velocity 
for the fixed X\ in a halo is the same irrespectively of the place of x±. Since we consider 
the strongly non-linear limit (x — > 0), we can neglect the effect of the boundary of the 
halo, which would be appeared when a?i is placed near the boundary. Then, we can change 
(■■■)iocai into the average for any fixed X\. Here we put the X\ on the center of the halo. 
Hence (...)i OC ai is exchanged into (...) x in which X\ is fixed to the center of the halo as we 
defined in eq. (|36D . Furthermore we find easily that (— • V\) x = 0. 

From eq. (|32|) , the density p is written by the following form, 



p oc a^ 5+ ")/("+ 3 )-V 



(38) 
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We put eq.(|38D into eq.(|33D, and then we obtain 



oti S f ream [n + 6) (oq\ 

Z -\- llstreamy 1 ' T <Jj 

It is found that when h stre am = (hi oca i = 0), the index of the density e must be 0. 
Furthermore, when h stre am = 1, the value of the index corresponds to that shown in Syer & 
White(1997). Syer & White(1997) shows the model that the density profile is generated by 
tidal stripping of small halos as they merge with larger objects. 

In our model, we can estimate h or the power index of the two-point spatial correlation 
function as a function of initial condition n. 

In the case of 3/2 < e < 3, as we can see in §3, we can consider that the two-point 
spatial correlation function obeys the power-law. Furthermore, when the self-similarity is 
satisfied, the power index of the two-point spatial correlation function can be expressed 
by eq.([L3|). On the other hand, when we assume that hi oca i is equal to 1, the index of the 
density profile e is expressed from eq.(^) by the following form; 

e = *i±2. (40) 

n + 5 

Here n must be greater than —1 because 3/2 < e < 3. On the other hand, from eqs. (|13D 
and flTTD, e is given by 

7 + 3 = 3h(n + 3) + 3 
2 2 + h(n + 3) [ ' 

From eqs. (|40"|) and fl4lf) , we obtain the form of the velocity parameter ftasa function of 
the index of the initial power spectrum n for n > — 1. 

h = -£±1 (42) 

2(n + 3) V ; 

We put eq.(pE2"D into eq.flTBD, and then we obtain 

7 = *££. (43) 

n + 5 

In the case of the < e < 3/2, the index of the two-point spatial correlation function 
in the non-linear limit (x — > 0) is equal to because £ = c\ + C2(x) 3_2e with 3 — 2e > 0. 
Furthermore, the velocity parameter h is equal to in the nonlinear limit (h = d,2X 3 ~ 2e ). 



-15- 



On the other hand, when h stream = h[ oca [ = 1 is assumed, the index of the density profile is 

e = 3(n + 3)/(n + 5). 
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